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Semiclassical expansions for traces involving Greens functions have two contributions, one from 
the periodic or recurrent orbits of the classical system and one from the phase space volume, i.e. the 
paths of infinitesimal length. Quantitative calculations require the control of both terms. Here, we 
discuss the contribution from paths of zero length with an emphasis on the application to Franck- 
Condon transitions. The expansion in the energy representation is asymptotic and a critical param- 
eter is identified. In the time domain, a series expansion of the logarithm of the propagator gives 
very good results. The expansions are illustrated for transitions onto a linear potential and onto a 
harmonic oscillator. 
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I. INTRODUCTION 

Many quantum properties, including the density of 
states and Franck-Condon transition matrix elements can 
be expressed as a trace of the Greens function times some 
operator Semiclassical expressions for such quanti- 

ties naturally divide up into two parts, one due to "paths 
of zero length" and one due to the longer recurrent orpe- 
riodic trajectories of the associated classical system 
The numerical and conceptual difficulties associated with 
the periodic orbit part have received considerable atten- 
tion in the literature (see the contributions to [f^g and 
the review Q). Zeta functions have helped to overcome 
many of these problems, at least in certain, well-behaved 
situations |9H12|| . More recently, higher order corrections 
to the dominant semiclassical contributions, in particu- 
lar in the neighborhood of caustics and bifurcations have 
been addressed |p^-fi8|. In applications to the photodis- 
sociation of molecules one has an additional source of 
corrections connected to the fact that the operator is a 
projection on the initial state and hence singular in the 
semiclassical limit. A way to deal with this has been 
proposed by Zobay and Alber Jll| . 

In mainly direct reactions the largest part of the cross 
section comes, however, from the paths of zero length. 
This leading order term is also known as the Thomas- 
Fermi contribution for the case of smooth systems or the 
Weyl-term for billiards. It measures the volume of the 
energy shell in units of Planck's constant, raised to the 
power of the number of degrees of freedom present |&Q . 
It turns out that in many situations one is too far away 
from this semiclassical limit (typically one needs a higher 
density of states) and so has to go beyond this leading 
order term. In billiards, the approximations to the den- 
sity of states regularly contain the subdominant contribu- 
tions from the surface and corner corrections, and even in 
smooth systems the leading order term alone will not do. 
The expansion in decreasing powers of energy or wave 
number, however, can typically be asymptotic at best. 
Building on their previous developments in the theory 
of asymptotic series, the behavior of the expansion has 



been illustrated by Berry and Howls |2C]] for the case of 
billiards: the series expansion indeed diverges and the 
rate of divergence is determined by short real or imagi- 
nary orbits of the classical system. 

In view of this it is perhaps not surprising that in 
semiclassical calculations of photodissociation cross sec- 
tions one also has to go beyond the leading order terms 
and that one encounters the same kind of divergences 
^|, pl|JTT| . For practical applications the problems then 
are how to estimate the importance of the higher order 
terms (without calculating them, of course) and how to 
improve on the series expansion. We will show here that 
comparison to a simpler problem, namely excitation onto 
a linear potential, suggests a useful parameter. For the 
second problem we analyze in some detail the behavior 
of three different approximations to the background term 
and identify the most useful one. For the sake of simplic- 
ity in notation all our calculations will be confined to one 
degree of freedom only. Generalizations to more degrees 
of freedom are straightforward. The main ideas will be 
illustrated for transitions onto a linear potential and onto 
a harmonic potential. Applications to photodissociation 
of water will be given elsewhere |3fJ . 

Various theoretical aspects of the classical and semi- 
classical limit of Franck-Condon transitions have been 
discussed previously in the literature. Much of the his- 
tory is reviewed in the paper by Dowling et al. [^2| , where 
also an interesting alternative phase space interpretation 
for Franck-Condon transitions can be found. Applica- 
tions to molecules can be found in the books p3|]24| ]. Of 
particular relevance for our discussion is a paper by Heller 
Q which contains background information as well as a 
discussion of the first few correction terms; we summa- 
rize some of his results in section IIA. Approximations 
in the time domain (which we take up in section 3) have 
been discussed e.g. in [p5|-p7[. 

The outline of the paper is as follows. In section 2 
we discuss the Wigner- and the Grammaticos and Voros- 
expansion in the energy representation. The behavior at 
large orders for a linear potential is analyzed in section 3. 
In section 4 we study approximations in the time domain, 
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which can then be connected to the energy domain by a 
Fourier transform (perhaps to be evaluated numerically) . 
The quality of the approximations is illustrated for the 
harmonic oscillator in section 5. Some conclusions will 
be drawn in section 5. 



It is interesting to note that under the trace the Wigner 
transform of the product of two operators maps into the 
product of the Wigner transforms. Generally, for a prod- 
uct of operators one has 

[AB] w (p,q) = [i] w (p,g)exp{|A} [B] w (p,q) (10) 



II. SEMICLASSICAL FRANCK-CONDON 

FACTORS with the differential operator 



A. The Wigner series 

We consider transitions from an initial state Vf^, typ- 
ically a Gaussian, to a manifold of final states at 
energy E. The system has one degree of freedom and 
the classical Hamiltonian on the upper potential energy 
surface is given by 



H(p,q) = 



2 m 



V(q). 



(D 



Quantum operators will be denoted with a so that e.g. 
the quantum Hamilton operator will be H. The Frank- 
Condon factors we want to calculate are the squares of 
the transition elements l0,E3|,E4f, 



p(E) = \(y i \H> E ) 



(2) 



Using |*_e)(*b| = S(E - H), where H is the Hamilton 
operator for the final electronic state, we can write the 
Franck-Condon factor as 



P {E) = {9i\8(E - £)\9t) 



(3) 

dx(*i\x)(x\8(E-H)\*i) (4) 

= trS(E-H)tl, (5) 

which has the above mentioned form as a trace over a 
Greens function, 



5(E -H) = lim G(E + ie) 

7r e^O 



— Im lim 



1 



e-0 e - H + ie 



(6) 

e > (7) 



times the projector II = |\&j)(\I>j| onto the initial state. 
Taking the Wigner transform of this expression one ar- 
rives at classical phase space traces over the Wigner 
transforms of the operators involved |2q] , 

i r+ca r+co 

p{E) = - / dp dq [5{E - H)] W (P, q) [fl]w(p, q) ■ 

(8) 

where the Wigner transform of an operator A is given by 

IS 



[A) w (p,q)= f dxe- ix *>l n (q + 



A 



q 



(9) 



A = 



v, 



v 9 - 



(11) 



where the differentiations act to the left or to the right as 
indicated. If the Hamilton operator is simply replaced by 
the classical Hamilton function, the leading order term 
for the density of states, the microcanonical measure on 
the energy shell, [S(E — H)]w ~ S(E — H), is obtained. 
The fact that formation of products and Wigner trans- 
forms does not commute is the source of quantum cor- 
rections which we want to study. 

A systematic procedure to calculating an expansion of 
the density of states in powers of h was suggested long 
ago by Wigner (2^] and applied to the Franck-Condon 
problem by Heller It uses the statistical operator 
P = cxp(— (3H) and its high temperature expansion for 
j3 — > 0. The statistical operator satisfies 



9/3 



HP = i ( HP 



PH 



(12) 



The symmetrized version on the right hand side is par- 
ticularly well suited for application of the Wigner trans- 
formation as it immediately shows that the quantum cor- 
rections come in even powers of Ti only. Using (10) for the 
product of two operators the equation for [P]w becomes 



- ±[P] W = [H] W cos {| A} [P] w . 



(13) 



The first few terms can be calculated by substitution of 
the ansatz 



[P] 



w 



-8H 



(l + h 2 c 2 + % a 



(14) 



This expansion contains only even powers of fi because of 
the symmetrization in (|lj) . Each coefficient c n is itself a 
polynomial in 0. For later reference we quote C2 for the 
Hamiltonian (111) MM. 



C2 = -Ff-I + /3 3 /3 



with 



h = —v" 

JZ 8m 
J 24m 



1 V"p 2 



24m 2 



(15) 

(16) 
(17) 
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where the primes denote derivatives with respect to po- 
sition. 

The statistical operator is the Laplace transform of the 
density of states, so to get back the quantity needed in 
the phase space trace, we need to do an inverse Laplace 
transform (see for more details). The leading expo- 
nential then maps into the delta function on the classical 
energy shell, 5(E — H), and the powers of j3 map into 
derivatives of the delta function with respect to energy. 
The resulting expansion thus becomes 



[S(E - H)] w = 5(E - H) 



^ 2 f2^S(E -H) + h'h^S(E - H) 



2 r d S 



■dE 



dE 3 



+h 4 



(18) 



Heller Q noted that instead of expanding in the form of 
an exponential times a power series in (3, one could al- 
ternatively expand as the exponential of a power series. 
For the terms given above this leads to an Airy function 
approximation, 

[6(E ~ H)] w = exp {-{H - E)f 2 /3f 3 - 2h 2 fl/27f 2 } 
■aAi(a(H -E + h 2 fl/f 3 )) (19) 



with 



a = (3ft 2 / 3 r 1/3 . 



(20) 



Up to and including terms of order h the delta- function 
expansion of the Airy function in (19) is equivalent to the 
Wigner expansion (|l8|). We will come back to this Airy 
function approximation below. 



B. The Grammaticos and Voros expansion 

The Wigner method of the previous section is a bit 
tedious when it comes to calculating higher oder terms. 
A very convenient algebraic method of expansion was 
developed by Grammaticos and Voros |30|. The Dirac 
measure [S(E— H)]w is expanded around the identity op- 
erator I times the classical Hamilton function, H (p, q) ■ I. 
The resulting series contains powers of the deviations 
H — HI, which have an explicit h dependence because of 
the Wigner equivalent of the product of two operators, 
Eq. [1^. To obtain this expansion, start from the integral 
representation, 



S(E - H) 
1 

2nh 
1 

2nh 



: iEt/h e -im/h dt 



,i{E-H)t/h e -i(H-HI))t/n dt 



(21) 
(22) 



_L / J(B-H)t/ty 



r=0 



(-ay 



(H-HI) r dt. (23) 



Next, take the Wigner transform on both sides and in- 
terchange summation and integration, 



[S(E - H)] w 

1 Y> (-i) r 

oo 

E 



(h - my 



w 



L-^6^(E~H(p,q))g r (p,q,h) 



r=0 



The derivatives of the delta function are with respect to 
the energy (as in (|l^)) and the universal coefficients Q r 
are the Wigner transforms of powers of the Hamiltonian, 



g r (p,q,h)= (H-H(p,q)IY 



w 



(26) 



They are universal in that other functions of the Hamil- 
tonian can be obtained by integration over energy, e.g. 



(27) 
(28) 



[f(H)] w = / dE f(E) 8{E - H) 

J — OO 

oo 1 

= Y J -f {r) {H{p,q))g r (p,q,h). 



and will contain the same coefficients. We will use this 
transformation for the statistical operator and the Greens 
function below. 

In principle, the Q r could be calculated by straightfor- 
ward expansions of the powers, and subsequent Wigncr- 
transformation. However, Grammaticos and Voros |3(]] 
give a more efficient algebraic technique, using Greens 
functions. They consider the transformation (|28|) with 
f{H) = l/(H — z) and thus find for Greens function 



(29) 



The Q r can be read off from an expansion of [G] w(j>i 9> K) 
in inverse powers of (H — z)~ r ~ l . The Greens function 
satisfies (H — z)G(z) = I. Using the symmetrized version 
as in (|l2|) and (|l3|), this becomes in the Wigner represen- 
tation 



[I] w = 1 = [(H- z)G(z)} 



ir 



[(H-z)] w co8[ - A) [G] w (z) 



= (H - z)[G] w + h 2n H 2 n[G] 



w 



(30) 



(31) 



The expansion in powers of Ti in the last line defines dif- 
ferential operators Tt n (p, q, V p , V,j). For n = 2, 4 and 
the Hamiltonian (PD) they are 



Ui = — I — 

= 



1 (}_&_ d 2 V{q) d 2 
8 \ to dq 2 dq 2 dp 2 
1 d 4 V(q) d 4 



384 dq 4 dp 4 



(32) 
(33) 



(24) 
(25) 
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The Wigner transform of Greens function can be ex- 
panded in powers of fi as well, 

oo 

[G] w (p, q,z,h)= ^2 h m G m (p, q, z) , (34) 

m=0 



and comparison of like powers of ti in (30) then gives 
recursion relations for the coefficients G m , 



Go(p,q,z) 



1 



H-z 1 
G2n{p,q,z) = - — 



— (^2^2kG 2 (n-k)^ 



(35) 
(36) 



Odd powers in h vanish because of symmetrization (as in 

(TO). 

Putting everything together one notes that by this 
method the Greens function becomes a double series, or- 
ganized in powers of H and in inverse powers of (H — z). 
Expression (|34|) emphasizes the former and (^9|) the lat- 
ter, so that the Q r can be obtained from ( |36| ) by collecting 
contributions with the same powers of (H — The 
structure of the expansion is such that there are always 
finitely many contributions to the coefficients Q r (for an 
example see the expansion for the harmonic oscillator in 
( |94| ) below). The first few coefficients are 



0o = l, 



Gi = 0, 



h 2 d 2 v( q ) 

4m dq 2 



(37) 



(38) 



(39) 
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h 2 f dv(g) y h 2 P 2 d 2 v(q) mr d 4 v( q ) 

4m \ dq J Am 2 dq 2 64m 2 dq 4 



This model can be solved analytically and can be used as 
an approximation for transitions in an arbitrary potential 
if a is taken to be the slope of the upper potential en- 
ergy surface at the maximum of the Gaussian initial state 
(which here is centered at qo = 0). Lengths q — q/X can 
conveniently be measured in units of 



A = (h 2 /2ma) 1 ^ 



(43) 



and energy in units of aA, i.e e = E/aX. The eigen- 
functions for this potential are Airy functions, and A 
sets the scale for the width of the oscillations near the 
turning point; the other oscillations in the wave function 
have shorter wave length. The Franck-Condon matrix 
elements can also be calculated exactly, 



P(e) = 



1 



e -fs 2 / 2 Ai(-(q + e))dq 



They depend on energy and on the ratio 
s = A/cr 



(44) 
(45) 

(46) 



of the length scale of the continuum wavcfunction to the 
width of the initial state. This parameter also contains h 
and thus indicates how 'semiclassical' one is. For small s, 
i.e. small h or large a and a broad initial state, the tran- 
sition integral will average over many of the oscillations 
of the Airy function and one can expect the leading clas- 
sical term to be reasonably accurate. However, for large 
s and narrow Gaussians, the initial state will probe every 
fine detail of the Airy function and the classical approx- 
imation will presumably not work well. We will come 
back to this point in section [VB. It is our aim now 



to show how these quantum expectations are reflected in 
\ 'the semiclassical expansion. As a first step we need to 
calculate the terms in the series (E5|). 



Given these expansions, we now investigate their large 
order behavior, first for the linear potential already stud- 
ied by Heller @. 



B. Asymptotic expansions for a linear potential 



III. LINEAR REPULSIVE POTENTIAL 

A. The quantum case 

The simplest model for a Franck-Condon transition 
onto a dissociating potential describes transitions from 
a Gaussian initial state 



1 _ e -(x~x ) 2 /2a 2 



onto a linear potential with Hamilton operator 



H = -—A -ax. 

2m 



(41) 



For the Hamiltonian ( f42| ) with its linear potential, 
equation ( |3l| ) contains only a single differential operator, 

H 2 = -\{H-z)K 2 
8m dq 2 ' 

all higher operators vanish. Then the recursion relations 
for the G2n reduce to 



Gin — —GoH.iG2(n-\) 



(47) 



( 42 ) Since 
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Go — 



H-z Pi - an - z 

2m ^ 



is linear in q, the action of Ti.2 on powers of Go is 
1 



(48) 



H 2 G% = -—a 2 n(n + l)G™ 
8m 



n+2 



(49) 



Together with the pre- factor —Go the coefficients G^n 
become 



_ / a 2 \" (3n)! 3 „ +1 
Comparing the expansion 



(50) 



[G]w = £ 



71=0 



a 2 \ n *2„(3n)! i 



8m / 3"n! (F - z) 



3/j 



+T (51) 



to the ?i - expansion of [G]vf in Eq. Q2S|) one can easily 
read off the coefficients Q r . The asymptotic series ( p5| ) 
for the Dirac phase space measure then becomes 



[8{E-H)] W = Y,- 



f /V?^" 



n! V 24m 



8 {Zn \E- H{x,p)). 



(52) 



This may be compared to the exact form of [8(E — 
H)]w computed from the exact propagator, 



[S(E H)] 



vv 



= Ai 



E - -kv 2 + aq \ I 



a 2 h 2 
8m 



\ 8m J 

1 f°° f 2 ft 2 1 

-y^^expj^-^-^^ (53) 



where z is a dimensionless integration variable. The se- 
quence of delta functions and their derivatives may be 
obtained Q by expanding the exponential with the z 3 
term in a power series and interchanging integration and 
summation, 

[S(E - H)) w 

i 00 i /n 2 h 2 \ n r°° 
1 ^1.(1^.) I dz(-iz 3 re^ E - H > 



2vr ^ n\ V 24m J 

n=0 v 7 



y__(_A) 5^(E-H) 
^ n\ V 24m / v ; 



n=0 



(54) 



These formulae can now be used for the calculation of 
the Franck-Condon factors. 



C. Asymptotic behavior 

The Wigner transform of the initial Gaussian is a 
Gaussian in phase space, 



= 2ex P 



(q - Qo) 2 



a 2 p 2 



(55) 



which helps to limit the domain of integration in (gj) espe- 
cially in cases of dissociation processes where the classical 
phase space is unbounded in certain directions. The ex- 
pansion for the density of states, eq. fl54|), then gives for 
the series expansion of the Franck-Condon factor 



n— 



(56) 



with the integrals 
I n {£, s) = 



\2 2 2/2 
-€) S —X f S 



H 3n ((x 2 -e)s)dx. (57) 



Just as the exact quantum result the expansion contains 
only a single parameter, s = X/a. The behavior for dif- 
ferent values of s at e = 0, i.e. for the maximum of the 
packet, is indicated in Fig. [l| For small s the terms de- 
cay rapidly and up to high n, but eventually they start to 
grow and to diverge rapidly: this is the behavior expected 
for an asymptotic series. As s increases, the turnover to 
divergence comes for smaller n and for s near 1 all higher 
order terms are larger than the first one. In Table [| we 
list the index of the smallest term as well as its size for 
several values of s. 

To estimate the rate of divergence of the series, one 
can substitute the generating function for the Hermite 
polynomials and express the integrals (p7h for e = as 



In(s) 



g3n 



dxe~ {x e) s x I s + 2t O £ ) s " f 



This integral cannot be solved exactly. In saddle point 
approximation for small s one finds just a single saddle 
(for some e there are three). The results is 



In{s) 



Q3n 

dt 3n 



2 VI - 2is 3 
Expanding in powers of t, one finally arrives at 

(-1)' 



(59) 



2l+r=3n 



11 



(2sT 



-1/2 
r 



(60) 



clearly showing a leading order behavior for small s of 
the form 



Us) 



(-1) 
(-1) 



3n/2 V£_(3n)!_ 
2 (3n/2)! 6 

(3rt+l)/2 V^ (3ra)! 

2 ((3n-l)/2)! 



n even, 



s n odd. 



(61) 



(58) 
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Together with the other prefactors in (|5^) a typical term 
in the sum for n even is 



(_l)8n/a ( 3n )| 
V7aA 12" n!(3n/2)!' 



„3n+l 



(62) 



and similarly for n odd. With the help of Stirlings for- 
mula, one can estimate the size of this term to be of 
the order [\/6(n/12 e) 1 / 3 s] 3n . The terms start to grow 
once the square bracket becomes larger than one, i.e. for 
n > n c m (2e/v6j s -3 - This strong s dependence is in 
accord with the variations in Fig. |l|. 

The information gained from this analysis of the linear 
potential can be used in more general settings: for excita- 
tions into some arbitrary potential V(q), the parameter 
a can be estimated from the slope of the potential at the 
maximum go of the initial state, a — —V'(qo). If the 
parameter s thus formed is small, the classical approxi- 
mation will be good (and the corrections small) , whereas 
some uniform approximation has to be tried if s is large. 



IV. TIME-DOMAIN APPROXIMATIONS 

A. The formal theory 

Another approach to the calculation of the mean den- 
sity of states exploits the relationship between energy- 
and time-domain. Since the mean parts are obtained by 
averaging over (large) energy intervals, they correspond 
to short time properties. This is the way the background 
is obtained in numerical calculations: the wave packet is 
propagated for a time long enough so that the overlap 
with the initial wave packet is minimal but short enough 
so that no major part has returned, yet. This finite time 
window is then Fourier transformed to obtain the back- 
ground term. 

This can be mimicked in the semiclassical approxima- 
tion. The Franck-Condon factors are related to Feyn- 
man's propagator in phase space by 



P (E) = -^-Re 
irh 



dt 



dqdp i 



Et/h 



[Tl]w(p, q) [K] w (p, q) . 

(63) 



where [a]vc i s the Wigner transform of the propagator 



[K]w(p, q) 



dxe- lxp/h (q + 



-mt/h 



(64) 



In the standard semiclassical approximation the smooth 
part of the Franck-Condon factor is obtained by use of 
the short time propagator from q — x/2 to q + x/2 for 
Eq. (|64|). In the simplest approximation (e.g. Berry and 
Mount 1972), K is replaced by a propagator for a free 
particle in a constant potential V(q), 



K { q + X ,2,q- X ,2,t) = (££) 



V /2 e K%X 2 -V(q)t)/h 



(65) 



The dominant part of the Franck-Condon factor as a 
function of the energy will come from that region in en- 
ergy where the turning point of the classical motion on 
the upper surface is near to the maximum of the initial 
state. The ultimate form of this is the reflection prin- 
ciple, where the wave function on the upper surface is 
replaced by a delta function so that the cross section is 
obtained by 'reflection' of the initial state on the poten- 
tial energy surface (see p|p3| for a discussion of its origin 
and its limitations). This implies that the straight line 
propagator on a constant potential is a rather poor ap- 
proximation since it does not account for a turning point 
nor the exponential damping when entering the forbidden 
region (See Fig. ||). 

To obtain a better approximation linearize the poten- 
tial around the center point q of the propagator, 

V(q') = V(q + {q' - q)) « V(q) + V(q) (q' - q) , (66) 
and use the exact propagator for the Lagrangian 
Ctf, «0 = E 2 - V '^) l' V ^) + V'(q) q , (67) 



viz., 

K (q + x/2, q-x/2,t) = 
(OT) 1/2 c*p(i[§z 2 -^)-^ 3 ]) ■ (68) 

The integration over x and time gives an Airy function 
phase space density, 



n,iE) I ^ &M(a(H(q,p)-E)) [IL]w(q,p) (69) 



with 



fi 2 |v"(g)l 2 N ~ J ; 

8m 



(70) 



Comparison with the Airy function approximation intro- 
duced by Heller, ( |l9| ) and (|2^), shows that the two coin- 
cide for a linear potential: then V" = as well as fi = 
and the two scales ( p0| ) and (|7(]) agree. 

The above approximative formula for the direct part of 
Franck-Condon factors is exact for a linear potential and 
contains an interesting limit case, the one of an almost 
flat potential. 
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B. Recurrence resonances 

If the potential vanishes, a — in (ft"2[), the Gaussian 
simply spreads and the Franck-Condon factor becomes 



pfree(E) 



2ma 2 e - 2mEa l % 



(71) 



When a = is substituted in (|69|), the scale a becomes 
infinite and the limit is singular. Similarly, the parameter 
s goes to infinity, indicating that the transition proceeds 
via the full quantum regime. Note that in the series ex- 
pansion ( p2| ) the transition can be performed easily, the 
coefficients of the derivatives of the delta function van- 
ish. This shows that the effect to be discussed here is 
non-perturbative. 

On the classical side, the limit is also singular in a 
certain sense. Fix the energy and consider the classical 
trajectories starting at some point qo- For vanishing po- 
tential, there are two straight line trajectories running off 
to positive or negative infinity, depending on momenta. 
As soon as there is a sloping potential, no matter how 
small a, the topology changes (see Fig. ||). For negative 
a, the trajectories running off to plus infinity will pick up 
in speed but be otherwise unchanged. The ones going to- 
wards negative infinity, however, run up to the potential 
and will have to return after some time. The position of 
this turning point, q c = —E/a if E is the kinetic energy 
at go j moves further out the smaller a and also the time 
for the trajectory to return increases with decreasing a. 
The consequences of this are that the return of the classi- 
cal trajectory introduces a long time scale, which, upon 
Fourier transform to the energy domain, will manifest 
itself in a modulation on small energy differences. 

These modulations are observed in the exact Franck- 
Condon factor (^||) in the case of large parameter s, as 
then the structure from the Ai 2 -part dominates (See Fig. 

A WBK quantization allows to connect the oscilla- 
tions with recurrent orbits. 

The positions of the 'resonances' in the exact Franck- 
Condon factor eq. (^) are given by the zeros of the first 
derivative of the Airy function which are approximately 
described by [ pl| : 



(?(• 



2/3 



i S - 4 
4 



(72) 



If we quantize semiclassically the action of the closed 
(recurrent) orbit which starts in the center of the Gaus- 
sian, is reflected by the potential and returns to the initial 
point with reversed momentum by requiring that 



(¥(• 



2/3 



(74) 



The quantization assumes orbits starting and ending in 
the center of the Gaussian and hence becomes more ac- 
curate for a narrow Gaussian. 

In the limit s — > these oscillations are suppressed: 
The initial Gaussian has a momentum distribution 
e~P s which is more concentrated around zero for small 
s. If we regard the initial wavepacket as a cloud of parti- 
cles with this momentum distribution, it will move quite 
uniformly for small s, whereas for large s the cloud will 
spread quickly. The fast moving parts are reflected and 
interfere with the slowly moving parts. As the particles 
obey e = p 2 — q, their energy distribution in the center of 
the Gaussian, q = 0, is e~ es . This energy distribution 
factor is responsible for the damping of the oscillations 
in the full cross section eq.(E3). 

If we take a look at the exact autocorrelation function, 



C(t) = 



1 



-t 2 / is 2 -it 3 / 12 



VI + its 2 

we have in the limit of the reflection principle 



C(t) 



-t 2 /4s 2 



(75) 



(76) 



which results in a Gaussian cross section, 
limit 



C(t) 



1 



Vl + its 2 



-it 3 in 



In the other 



(77) 



the first part resembles the autocorrelation function of a 
freely moving spreading wavepacket, but the second oscil- 
lating term accounts for the recurrent path contributions 
to the propagator. 

It is therefore reasonable to term the oscillating struc- 
tures 'recurrence resonances'. Although they are not con- 
nected to periodic orbits which imply infinitely many 
revivals they are associated with recurrent orbits, and 
their autocorrelation function differs from that of a sim- 
ply spreading wavepacket. 

Already Condon J3^] and then later Hunt and Child 
p3[ and Tellinghuisen |34| have mentioned these oscil- 
lations to appear on very shallow dissociative potentials 
and described them with the help of an opposite reflec- 
tion principle, namely that the continuum wave function 
is reflected in the cross section |3Jj : For large s, we have 
in the Franck-Condon factor eq.^44): the initial state cen- 
tered at q = which acts like a needle scanning the con- 
tinuum state: 



S(E } - - 



-sX 



^2m(E + ax)dx = 2n{n + ~) (73) 



(note that there is just one turning point, hence the con- 
tribution 1/4 instead of 1/2 from the Maslov indices) we 
obtain resonance positions 



2vrcr 



and the cross section becomes 

p( 6 ) = ^Ai'(-*) 



(78) 



(79) 
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This example nicely illustrates the importance of recur- 
rent orbits since they are responsible for ALL the struc- 
tures which often would be associated with periodic or- 
bits. 



C. Higher order approximations in the time domain 

In both the energy and the time domain, the Airy func- 
tions are intuitively appealing approximations based on 
the first few correction terms, but it is not obvious how 
to systematically improve on them. A most convenient 
method uses the Wigner expansion of the time prop- 
agator, obtained from (^8|) with the Grammaticos and 
Voros method and the universal coefficients Q r , compare 
eq. (p5|). The Wigner transform of the propagator has 
the expansion 



,-iHt/H 



oo 
r=0 



iHt ' h g r {q,P,h) , (80) 



with the same Q r as calculated before. This means that 
the autocorrelation function of the wave function, 



C(t) = (*i|¥ r ) = tr e- iAt/R fl 
dpdq 



-iHt/H 



W 



0) q) [n]w(p, q) 



(81) 
(82) 



can be constructed with the help of this series. The p- 
integrations can always be done analytically and we are 
left only with the coordinate integrations. The cross sec- 
tion is then obtained by a Fourier transform, 



9 f°° 
p(E) = -Re J e iEt / n C(t) X (t), 



(83) 



where x(t) is a window function that confines the region 
of integration to the short times needed for the direct 
part. This window is moreover important as the terms 
in the expansion ( |S0| ) may grow without bound for large 
time, influencing the convergence of the integrals. As 
soon as terms beyond the linear potential approximation 
are included the autocorrelation function may show re- 
currences structures which have to be switched off by 
means of a window. 

Eq. 82 has the advantage over ( |l8| ) that it avoids the 
energy derivatives which impose the highest demands on 
the numerical accuracy of the phase space integrations. 
The Fourier transform to energy can be done effectively 
using fast Fourier transforms. However, if the series ex- 
pansion ( |8C| ) is used directly, the high powers for large 
times have a devastating effect numerically and nothing 
is gained. A way out of this dilemma is suggested by the 
manipulations that lead to (p4[). Perhaps it is possible 
to sum the power series expression in t in (80) into some 
exponentiated form, the first term of which would be an 
Airy function. This may be achieved using 



ay 



Q r = exp 



-k 



k 



(84) 



The expansion coefficients Q r are known from the alge- 
bra in the energy representation and the coefficients J- k 
can be obtained from the Plemelj-Smithies recursion re- 
lations 



J 71 



n-n+l 



n-1 



n 



fe=i 



(n-k)\ 



Since for t — wc must have 



-i(H-H)t/h 



coefficient J-q = 0. The next few are 







T-2 








?4 





(85) 
1, the 

(86) 
(87) 



(89) 



The above formulas have an interesting phase space in- 
terpretation: According to Berry fjsfl we have now taken 
into account the fringes in phase space surrounding the 
energy shell. The difference to Berry's fringes-formula is 
that he omits the second term in t 2 on the assumption 
that the stability of the short trajectory segments does 
not change or that it changes slowly compared to the 
other terms linear or with time. He assumes the term 
quadratic in t, which is connected with the stability of 
the short trajectory segment, to be constant. 



V. MATRIX ELEMENTS BETWEEN HARMONIC 
OSCILLATOR EIGENSTATES 

To illustrate the quality of the various approximations 
we have to go beyond the linear potential (since there 
the Airy functions are exact) and turn to excitations into 
harmonic potentials. Contrary to expectations based on 
the usual close relationship between classical and quan- 
tum dynamics in harmonic oscillators, the semiclasscial 
expressions for the Franck-Condon transitions are not ex- 
act and thus a useful test of our formulae. 

As in the linear case, the initial state is a Gaus- 
sian (ffl] ) of width a and a Gaussian Wigner transform 
in phase space, (p5|) . The Hamiltonian now is 



H 



fc2 2 

n . mui 



2m 



-A 



(90) 



The length scale I — y/K/mui is characteristic of the wave 
length of the ground state, all higher excited states oscil- 
late on shorter scales. It thus plays the same role as the A 
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in the case of the linear potential, setting the largest scale 
for quantum oscillations. The Franck-Condon transitions 
depend on the energy and on the single parameter 



r = ija 



(91) 



which measures the size of the quantum oscillations rel- 
ative to the width of the initial state. 

The Wigner transform of the initial state has a maxi- 
mum in phase space near x — xq and p — 0. The leading 
order classical phase space average then yields a cross 
section which increases up to E max w murx\j1 and de- 
creases for higher energy. Since the quantum spectrum 
is discrete, the quantum transition strengths are discrete 
as well. A direct comparison between the leading order 
classical Franck-Condon spectrum 



Po(E) 



dqdp 
2nh 



S(E - lf(p J «))[ J |* < )(*<|] w (p, 9 ) (92) 



and the quantum data is shown in Fig. ^[ 

The quantum calculations show that for large h and 
large r the quantum cross section is shifted towards 
higher energies. 

The Grammaticos and Voros scries (^5|) for the har- 
monic oscillator can be calculated as before for the linear 
potential. The asymptotic expansion of the phase space 
trace can be carried out quite easily for the scaled Hamil- 
tonian H — fiu (|p 2 + ^q 2 ) where q — q/£ and p — pl/%. 
Again, there is only one operator in the series (j3l|), 



n 2 = 



1 d 2 



8 



)r,2 



8dp 



(93) 



and the others vanish 
then find up to order 0(h 



From the recursion relations we 
with the help of MAPLE: 



S(E-H)] W =5{E-H) 
- (he) 2 (~IS"(E -H) + ±E5"'(E H)^ 



v ; V 384 160 



■E 



6 (5) {E-H) 



46080 v 7 

sf jn_ 5 W(E-H) + -¥?-E8<n(E-H) 
16080 v ' 161280 v ' 

-E 2 S^(E -H) l ^—EH^(E - H) 

!0 v ; 11612160 v ' , 



v ; V 46080 v ' ' 161280 

- 11648 E^(E - H) - 1120 EW\E - H)) 
10321920 v ; 11612160 v 'J 

(Hto) 8 ( 1385 8&(E -H)- 4757 E5W(E - H) 
K ' V 10321920 v ; 11612160 v ; 

+ 1529 E 2 6^(E -H)- -^—EH^{E - H) 
6451200 v ; 414720 v ' 

+ — - — E 4 5^ (E ~ H)] (94) 
497664 v 7 



which involves derivatives of the delta function up to 12th 
order. Incidentally, this expansion shows rather clearly 
the structure of the two different series for Greek's func- 
tion, one in powers of Ti ( |34| ) and one in derivatives with 
respect to energy, ( p5| ) . Fig. || shows the quantum correc- 
tions for the harmonic oscillator Franck-Condon factors 
for different values of the expansion parameters. Two 
trends can be recognized: For large r, divergence sets in 
very early, even from the first term on, and the smaller 
the energy (or the smaller the local gradient of the po- 
tential) the worse the asymptotic expansion. In both 
cases the leading order term gives the correct qualitative 
shape of the Franck-Condon factor. For large energies, 
the higher order corrections are small and the series is 
close to the exact values. For lower energies, the higher 
order corrections increase and alternate in sign - the typ- 
ical features of an asymptotic series. The oscillations 
increase rapidly with r so that the series becomes prac- 
tically useless without re summation. 

The Airy correction (|6^) has been computed and gives 
remarkable results in the extreme quantum regime: Fig. |^ 
shows the lowest order classical approximation to the har- 
monic oscillator matrix elements and the Airy function 
approximation. This confirms the importance of turning 
point corrections and tunnel trajectories in the Franck- 
Condon region. 



VI. FINAL REMARKS 

We set out to calculate quantum corrections for 
Franck-Condon factors in the classical phase space trace 
(||). In a nutshell, there are three main results: i) For the 
expansion based on ( |25| ) in the energy domain the situa- 
tion seems to be that either the corrections are small to 
begin with (in which case one would be happy with the 
leading order result) or they are large and a resummation 
of the series is required as it is asymptotic at best, ii) 
In the time domain the expansion in the exponent, (|85|), 
provides a more useful representation and contains non- 
perturbative effects (the recurrence resonances) . iii) The 
quality of the approximation in all cases is controlled by 
a single parameter, which may be estimated from the ra- 
tio of a de Broglie wavelength to the width of the initial 
state as in @) and (|l]). 

Clearly, the parameter controlling the convergence can 
be put to immediate use and can help to identify whether 
a classical approximation will suffice. An example will be 
given in the photodissociation for water [ pdpS| ] . 

The relationship between the series expansion in en- 
ergy with its derivatives of delta functions, (^5|) and the 
exponentiated one in the time domain, (|S5|), is subtle and 
at present not fully understood. The numerical observa- 
tion is that it is in very good agreement with the quan- 
tum results and that no divergences seem to occur. The 
calculation leading to ( |54| ) for the expansion of the Airy 
function shows that at least in this case all divergences 
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are removed instantly: the exponentiated series stops af- 
ter the cubic in time. In general circumstances, this will 
not be the case and the series will continue. Then the 
question is whether the power series in time, (pQ), or its 
exponentiated version is convergent or whether at some 
late time, perhaps related to recurrent or periodic orbits, 
deviations between quantum and semiclassical time evo- 
lution become noticeable. For instance, in the case of the 
harmonic oscillator the propagator has an exact recur- 
rence after half a period which can only be accounted for 
by a divergence of the series. The success of the present 
calculation in the time domain would then be related to 
the fact that the time evolution is followed for a short 
time only. However, from a practical point of view, that 
is all one is asking for in the calculation of the direct part 
in Franck-Condon transitions, and we therefore propose 
to use ( |S2|) with the exponentiated expansion ( |85|) and a 
fast Fourier transform to the energy domain. 

The formulae presented allow for an accurate semiclas- 
sical calculation of the background term and thus for the 
largest contribution to the cross section in predominantly 
direct reactions. The positions and widths of resonances 
can be calculated from periodic orbit expansions and zeta 
functions. Calculation of the cross section however again 
requires improved formulae which take into account the 
neighborhood of the trajectories and the fringes in phase 
space are needed. This would then allow to calculate 
all aspects of cross sections semiclassically. Work along 
those lines is in progress. 
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FIG. 1. Series expansion for Franck-Condon transitions 
onto a linear potential for a Gaussian initial state. Shown 
is the decadic logarithm of the relative error of the series in- 
cluding the first n terms taken at the energy corresponding 
to the center of the Gaussian. 

FIG. 2. Contributions to the Wigner propagator in posi- 
tion space. The simplest approximation to the short time 
Wigner propagator uses the free one on a constant potential 
energy V(Q). This is appropriate for high energies like E\. 
Tunnel effects are important at lower energy E2 and can be 
included by use of the propagator on the linearized potential. 

FIG. 3. Contributions to the Wigner propagator in phase 
space. Two paths contribute to the semiclassical propagator 
from Q — X/2 to Q — X/2: the direct one from phase space 
point 2 to 3 and the path through the turning point from 1 
to 3. 

FIG. 4. Exact Franck-Condon factors for the transition 
onto a linear potential. For small s the cross section eq. ( fil] ) 
has a single maximum. For larger s recurrence resonances de- 
velop. In rescaled units they have a spacing given by eq. (fri)). 
In original units the oscillations become denser and denser. 

FIG. 5. Comparison between the n-th harmonic oscil- 
lator eigenstate matrix elements (n| 1 i'i)(^ r i|n) of the projec- 
tor onto a Gaussian ^i{q) and the correspondin g cl assical 
phase space average of the Wigner transform eq. (|92[). The 
center of the initial Gaussian remains the same for all fig- 
ures. The oscillator eigenvalues are E n — Tiuj(n + 1/2) with 
h = 4.0, 2.0, 1.0, 0.25 and w = 1.0. 



FIG. 6. Quantum corrections to Frank-Condon factors for 
transitions to harmonic oscillator eigenstates. The different 
curves include all terms of eq. (^) up to the indicated order 
in h. For larger energies, the corrections improve the lead- 
ing order result. For low energy they develop oscillations and 
diverge right away. 

FIG. 7. Comparison between the Weyl term, the uniform 
approximation and the exact results for transitions in har- 
monic oscillators. The exact result is shown by full circles, 
the Weyl term eq. as a solid line and the uniform ap- 

proximation eq. (p9f) as a dashed line. 

TABLE I. Asymptotics for the cross section of a Gaussian 
on a linear potential: We calculate the series up to its smallest 
term. The number of terms included grows with decreasing 
expansion parameter s. log 10 Ap(e = 0) is the logarithm of 
the relative error of the series including the first n terms taken 
at the energy corresponding to the center of the Gaussian. For 
large s > 1 the asymptotic series is practically of no use. 



s 


n 


log 10 Ap(e = 0) 


1.8 





0.03 


1.6 


1 


-0.16 


1.4 


1 


-0.59 


1.2 


1 


-1.38 


1.0 


1 


-1.68 


0.8 


5 


-2.50 


0.63 


21 


-4.51 
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